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Abstract 

This paper focuses on studying the fundamental performance limits and linear 
dispersion code design for the MIMO-ARQ slow fading channel. Optimal average rate 
of well-known HARQ protocols is analyzed. The optimal design of space-time coding 
for the MIMO-ARQ channel is discussed. Information-theoretic measures are used 
to optimize the rate assignment and derive the optimum design criterion, which is 
then used to evaluate the optimality of existing space-time codes. A different design 
criterion, which is obtained from the error probability analysis of space-time coded 
MIMO-HARQ, is presented. Examples are studied to reveal the gain of ARQ feedback 
in space-time coded MIMO systems. 

Index Terms: Space-time code, Hybrid ARQ, Multiple-antenna channels. Linear Dispersion 
Code. 



1 Introduction 



Hybrid ARQ (HARQ) techniques combine the automatic-repeat-request (ARQ) feedback 
with the forward error correction (FEC) codes to achieve better reliability and higher 
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throughput [1]. They are widely used in most of the contemporary communication sys- 
tems, as the delay tolerance of many data services allows using retransmission to recover 
erroneous packets. A good summary of the progress of HARQ schemes is presented in [2]. 
From an information-theoretic perspective, HARQ systems can be viewed as channel with 
sequential feedback, and there has been some research on the fundamental limits of this chan- 
nel. An information-theoretic throughput and delay analysis of several HARQ schemes in 
a Gaussian collision channel is presented in [3]. Throughput analysis of incremental redun- 
dancy HARQ in the block-fading additive white Gaussian noise (AWGN) channel is carried 
out in [4]. Optimal average rate performance of a scalar slow fading channel with ARQ 
feedback is done in [5,6]. 

Most of the aforementioned research is restricted to scalar channels. A natural question 
is how to incorporate the design of HARQ into multiple-input-multiple-output (MIMO) 
systems, and it has sparked many research activities. Most existing literature regarding 
MIMO-HARQ design falls into the following two categories. One is the information-theoretic 
analysis, for which the focus is on extending the Diversity-Multiplexing Tradeoff (DMT) 
framework [7] to consider ARQ feedback [8-10]. The other is the joint design of MIMO 
transmission and ARQ feedback. Receiver processing of ARQ retransmissions is studied 
in [11,12]. It has been shown that multiple ARQ packets combining is non-trivial in a MIMO 
setting, especially when soft-output is desired [13]. For the transmitter design, two different 
approaches have been taken to exploit the additional spatial degrees of freedom. One is to 
study the bits/symbols rearrangement for retransmissions to exploit diversity [14-17], and 
the other is to explore linear precoder design [18,19]. 

Despite the many efforts in studying MIMO-HARQ design, there are still some important 
problems left unanswered. From the perspective of information-theoretic analysis, the DMT- 
based approach focuses only on the high Signal-to-Noise Ratio (SNR) asymptotics and gives 
the tradeoff between multiplexing gain (pre-log), diversity gain (error exponent) and ARQ 
delay. It also requires a family of space-time codes (STC) whose rates scale with the SNR. 
In practice, however, one is also interested in designing a fixed-rate STC that operates well 
within a range of finite SNRs. It is unclear what is the optimal MIMO-HARQ throughput 
in this scenario, and how to design the optimal STC to achieve it. For the practical HARQ 
retransmission protocol design, the aforementioned works provide separate and ad hoc de- 
signs. There is no general framework that unifies the existing results and guides the optimal 
design. 

This paper is devoted to solving these problems. We begin by introducing the system 
model in Sec. [2l The performance limits of existing HARQ protocols in a MIMO setting 
are investigated in Sec. [31 Sec. H] studies the fundamental average rate limit of STC-based 
MIMO-HARQ protocols, and gives the optimal STC design criterion. Many existing STCs 
are re-visited under this framework, both analytically and numerically. A different approach 
based on the error probability analysis is given in Sec. 0, where the so-called n-th pairwise 
error probability plays a key role in the analysis. Optimal STC design criterion based on 
the error probability analysis is proposed and exemplary STCs are studied. Finally, Sec. [6] 
concludes the paper. 
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Throughout the paper the following notations will be used. Matrices and vectors are denoted 
with bold capital and lowercase letters, respectively. A{i,j) is the (i,j)-th element of the 
matrix A, and a(i) is the i-th element of the vector a. E[-] denotes the expectation of a 
random variable, a* denotes the conjugate of the complex number a, Re {a} is the real part 
of a, and a* is the transpose of vector a. We use Tr (A) and HAHi? to denote the trace 
and Frobenius norm of matrix A, respectively. A^ and conj (A) are the Hermitian and the 
conjugate of the complex matrix A, respectively. The Kronecker product operation of two 
matrices is denoted by ®, and vec (■) is the vectorization operation. Finally, |VV| denotes 
the cardinality of the set W. 



2 System model 

In this section, the MIMO channel model with ARQ feedback is described, and existing 
HARQ protocols that will be analyzed in Sec. [3] are introduced. 



2.1 MIMO channel with ARQ feedback 

This paper studies a single-user multiple-antenna wireless system equipped with Lt transmit 
and Lr receive antennas. It is assumed that there is an error-free and delay-free ARQ 
feedback link, which indicates to the transmitter the successful decoding by acknowledgement 
(ACK) and failed decoding by negative acknowledgement (NACK). The maximum number 
of ARQ transmissions is denoted by N, which is also called the maximum allowable ARQ 
rounds in the literature. With this constraint, if after a total of N transmissions the receiver 
still cannot decode the message successfully, no further attempt will be tried and a decoding 
failure is declared. 

It is assumed there is a very large pool of information messages available at the trans- 
mitter. As soon as the transmission of current message (including possible retransmis- 
sions) is done, the next message is encoded and transmitted immediately. Consider a set 
of uniformly distributed messages W. For the transmission of each message, the infor- 
mation message w E W is sent to a MIMO-HARQ encoder, which generates N matrix 
sub-codewords {X.n{w)}^^^, each corresponds to one ARQ round. Sub-codeword Xn(iy) will 
be transmitted in the n-th ARQ round if the decoding fails at the {n — l)-th round. Denote 
X("')(ti;) = [Xi(w), ■ ■ ■ , X„(w)], Vra = 1, ■ ■ ■ , A, which is the overall codeword after the n-th. 
ARQ round. The receiver performs optimum decoding based on all the received packets in 
ARQ round 1, ■ ■ ■ ,n. Each sub-codeword satisfies X„(tf) G C^*^'^"/^", where L„/r„ denotes 
the length of the sub-codeword, r„ is determined by the "coding rate" of the space-time 
structure as r„ = where Kn and T„ are the number of symbols and the time slots per 
space-time codeword, respectively, and L„ is the capacity-achieving channel code length. 



3 



With this model, the overall transmission rate in the n-th ARQ round is 



^(n) ^ l0g2 |W| 



Vn = l,---,iV. (1) 



The MIMO channel H G C^''^^* is assumed to be a random matrix with independent and 
identically distributed (i.i.d.) complex circularly symmetric Gaussian entries with unit vari- 
ance. The channel matrix H is assumed to be constant within the channel coherence time 
(Tc symbols), and changes independently to a different value in the next coherence block. 
Depending on the relationship between Tc and {L„/r„}^^^, different MIMO-HARQ mod- 
els have been studied in the literature [8,10]. In this work, a slow fading MIMO channel 
model is considered, where ^ '^n=i^n/i^n- This means the MIMO channel remains 
unchanged throughout all possible retransmissions. This models the worst-case scenario 
compared to [8,10], as there is no time diversity to exploit by retransmissions. 

With the described model, the received signal corresponding to the n-th transmission can 
be written as 

Y„ = y^HX„, + Z„ (2) 
and the overall received signal after the n-th ARQ round is 



YW = y^HX(") + Z("), (3) 

where Y*^") = [Yi, ■ ■ ■ , Y„] and X*^"\ Z*^"^ are similarly defined; the additive noise Z*^"^ has 
i.i.d. entries 2i^"'\i,j) ~ CA/'(0, 1). The overall transmit codewords ^^^\w) are normalized 
to satisfy the average power constraint^ 

\ , E||XWH|||<L, (4) 

where the expectation is over all possible codewords. Therefore, the average SNR per receive 
antenna is SNR. The receiver is assumed to have perfect knowledge of H. This is a reasonable 
assumption since the slowly varying nature of the channel helps receiver channel estimation. 
The transmitter is assumed to have no knowledge of H prior to the transmission, but it 
should be noted that ARQ retransmissions essentially inform partial CSI to the transmitter 
in a sequential fashion [6] . Fig. [1] illustrates the system model under consideration. 



2.2 Hybrid ARQ strategies 

HARQ incorporates FEC coding into ARQ feedback to enhance the retransmission perfor- 
mance. Depending on transmitter's retransmission strategy and receiver's decoding from 

^The transmit power constraint will be discussed in more detail when the optimal LDC design for MIMO- 
HARQ is presented in Sec. [H 
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Figure 1: The MIMO-HARQ system model. 

multiple packets, several efficient protocols have been proposed and studied. In this work, 
two traditional stop-and-wait HARQ protocols are considered: 

1) Incremental Redundancy (IR). This is a code combining scheme [20]. At the trans- 
mitter, information message is encoded into an overall codeword of length Yln=i ^n, 
and then the codeword is serially punctured into sub-codewords each with length Ln- 
These sub-codewords are mapped into the transmit matrix codewords [Xi, • ■ ■ ,Xjv]. 
At the n-th transmission, X„ is sent through the channel as additional (when n > 1) 
redundancy symbols to reduce the transmission rate. The receiver tries to decode the 
message based on all the packets it receives up until round n. 

2) Chase Combining (CC). This belongs to the category of diversity combining [20]. 
Upon each retransmission request, the transmitter simply repeats the same packet: 
X„ = Xi, "^n = 2, ■ ■ ■ , N. Different diversity combining schemes can be used at the 
receiver, among which the Chase combining [21], which essentially is a maximum ratio 
combining of all the received packets in the scalar channel, gives the best performance. 

These well-known HARQ protocols will be analyzed in a MIMO setting and the optimal 
performance is reported in Section [31 

Traditional HARQ protocols are designed to extract time diversity in a scalar fading channel. 
The performance depends critically on the channel variation over retransmissions. The 
slow fading channel model considered in this work allows no time diversity for the HARQ 
protocols to exploit, and hence packet-combining-based schemes have limited performance 
[6]. Multiple antenna system provides additional dimensions for spatial diversity, which is 
not effectively exploited by traditional protocols. On the other hand, although provides 
the best performance, IR has higher complexity due to the generation of rate-compatible 
channel codes or the rateless codes. Both problems call for the design of HARQ protocols 
that jointly consider the space-time structure to exploit spatial diversity and low-complexity 
ARQ retransmissions. This is the focus of Section H] and [51 

The figure of merit in this paper is the long-term throughput, and an information-theoretic 
view is taken for the MIMO-HARQ design problem. The throughput of a HARQ protocol is 
defined [20] as the average number of bits accepted by the receiver in the time it takes to send 
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the data packet. As discussed in [6], how to compute the throughput depends on the specific 
apphcations and the assumptions on how to use the channel. Typical derivation involves 
the evaluation of average transmission time. Renewal theory has been successfully applied 
to this problem [3,22]. However, this formulation leads to the problem that the empirical 
channel distribution does not match the true channel statistics, and the reason is that the 
length of channel uses is determined by the instantaneous channel state. This problem 
has been discussed in [6,23], and it has been shown [6] that for some applications, average 
rate [24-26] gives the long-term successful communication rate and is a more reasonable 
metric. The average rate is chosen as the performance measure in this work. 



3 Performance limits of existing HARQ protocols 

This section studies the optimal average rate performance where the existing HARQ proto- 
cols are directly applied in a vector channel. We will first derive the average rate for general 
HARQ systems, and then apply the result to IR and CC. The general average rate expression 
is a fundamental result and will also be useful in developing optimal space-time HARQ in 
Section HI 

Before analyzing the average rate, some important events are defined and the associated 
probabilities are studied. Due to the randomness of the channel matrix H and the limit N 
on the ARQ deadline, the successful communication rate i? is a random variable. Define two 
events for ARQ round n, Vn = 1, ■ ■ ■ , A^: 

An = Decoding at the end of ARQ round n is successful; (5) 
Sn = ARQ round n is activated. (6) 

There are two important observations about these events, both originated from an informa- 
tion theoretic viewpoint!^. First, a decoding failure is equivalent to a channel outage, i.e., the 
instantaneous channel capacity cannot support the current transmission rate. This is be- 
cause long-blocklength capacity-achieving channel coding "averages out" the additive noise. 
Thus one has 

Pr {An} = Pr (Hg)) > } . (7) 

where C''"^(Heq^) is the equivalent channel capacity at ARQ round n; i?*^"^ is the overall 
communication rate at round n, which is given in Equation ([1]). The second observation is 
that additional transmissions can only increase the successful decoding probability, via either 
increasing C''"^(Heq^) or reducing or both. This observation suggests 

An-l C An, (8) 

and hence 

Pr{5„} = Pt{A'i,-- - 

= Pr{AlT}. (9) 
^As one will see in Sec. O finite coding length leads to some different conclusions for the HARQ design. 
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These observations play an important role in analyzing the information-theoretic perfor- 
mance limits of MIMO-HARQ. 

With the definitions ([5]) and ([6]), the successful communication rate R can be expressed as 

/ if {5„,A},Vn = l,-- - ,iV 

"^=10, if:4^ 

if {XIT,X},Vn= I,--- ,iV 
0, if^liv ^ ' 

where the second equality comes from ([9]), and the convention Pr{^o} = is used. The 
average rate is then 

N 



R = ^i?Wpr{X-i,X} (12) 

n=l 
N 

= ^i?('^)Pr {C("-^) (H^:;-^)) < C(") (H(^)) > , (13) 



n=l 



where we define i?^"^ = oo (since the denominator in ([T]) is zero) and C^^-* (Heq) = 0. A 
different expression can be obtained if one notices that 

Pr [AT^, An] = Pr {An} - Pr { A-i} , (14) 

which is true since An-i ^ An- The average rate expression (fT2|) then can be written as 

N 



R = 5^i?(")(Pr{X}-Pr{X-i}) 

n=l 
iV 

= J] (i?^"^ - ^^"'"'^) Pr { A,} (15) 

n=l 
AT 

= - Pr (H(^)) > , (16) 



n=l 



where we define R^^^^^ = (since the retransmission stops after the n-th round). It should 
be noted that the term i?'-"-' ~ R^^^^^ denotes the rate decrease in the {n+ l)-th transmission. 

The equivalent channel capacity C^"-'(Heq'') at ARQ round n is determined not only by the 
instantaneous channel matrix H, but also by the specific HARQ protocol and the transmit 
covariance matrices {Kj}"^-^, where Kj is the transmit covariance matrix at ARQ round i 
and could be different over retransmissions. As one shall see in the sequel, these degrees of 
freedom shall be carefully exploited in the MIMO-HARQ protocol design. 

With the general average rate expressions (fT3|) and (fT6|) . we proceed to study IR and CC, 
respectively. 
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3.1 Incremental Redundancy 



The average rate for general IR can be obtained directly from (11611 : 

N 



R.„ = J2 (i?^"^ - Pr {Ci") (hI;),) > i?W } , (17) 



n=l 



and the remaining problem is to analyze the capacity of IR after each ARQ round. The 
capacity for a given MIMO channel H and total transmit covariance matrix K is known to 
be [27] 

C^i^o (H, K) = logdet (II, + HKH^) (18) 
where K has to satisfy the power constrainllfl 

Tr(K)<SNR. (19) 

Assuming transmit covariance matrix K„ is used in ARQ round n, the overall IR capacity 
where different portions have different covariance matrices is a TDMA-type one: 



X]j=l C*mimo (H, Kj) Li/Vi 



D(n) 



1=1 ^ 

where ( 12U]) is obtained from ([T]). 

The general IR average rate maximization problem can now be formulated as 

maximize Zti (^^"^ - R^^^'^) Pr {E7=i (w " W^) ^mimo (H, K,) > i?^} 
subject to > (21) 

Tr(K„) < SNR,Vn = I,-- - , 

Solving this problem is extremely difficult, which comes mainly from the search of optimal 
covariance matrices {K^}^^-,^. Evaluating the objective function in Problem 0211) requires the 
distribution of MIMO capacity Cmimo (H, K) as a function of the covariance matrix K, which 
is an unsolved problem. In fact, even the simpler problem of finding the optimal covariance 
matrix that minimizes the outage probability is still open. 

To make progress, this paper will focus on the isotropic Gaussian input distribution over all 
transmit antennas: 

K„ = K = ^1^,. (22) 



■^The constraint (fT9| assumes that each ARQ round consumes the same average power. One can also 
dynamicaUy aUocate power among ARQ rounds, which leads to a different power constraint and average 
rate expression. The case for scalar fading channel has been addressed in [6] . 
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This is by no means the optimal solution but very easy to deal with and widely assumed in 
literature. With fl22l) the IR capacity is the mutual information of a MIMO channel with 
isotropic Gaussian input distribution: 

= logdet(^Ii, + ^HH^^, (23) 
and the IR average rate maximization problem becomes 

maximize E1i (^^"^ - ^^"+'^) Pr {c.imo (h, ^^I^,) > R'^^A 
subject to > 



(24) 



For the remaining of this paper, the isotropic Gaussian input distribution is assumed unless 
the covariance matrix is explicitly specified. 



3.2 Chase Combining 

In CC, each retransmission effectively increases the number of receive antennas by Lj. via the 
repetition of previous packet. Thus after ARQ round n, the equivalent system is a {Lt,nLr) 
MIMO with an equivalent channel matrix 

H(^),, = 1„®H, (25) 

and 

C^:^ (aScc) = ^C7^imo (H^?,,, K) . (26) 

Meanwhile, since each retransmission is a simple packet repetition, the effective rate after 
ARQ round n is 

= -, (27) 

n 

where R is the rate in ARQ round 1. Hence, the general CC average rate maximization 
problem is formulated as 

maximize Zn=i f {Cmimo (h^'c^ k) < i?, C^i^o (h^cc, k) > i?} 
subject to R>0 (28) 
Tr (K) < SNR. 

Note that due to the simple packet repetition in retransmissions, only one covariance ma- 
trix K and one rate parameter R need to be optimized. If the isotropic Gaussian input 
distribution is assumed, the problem becomes (-R > can be dropped): 

maximize f P^ (h^-c^ ^I^*) < ^' (^5?-' ^I^*) ^ ^} " ^^9) 
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Further simplification is possible if one notices that 

(hSJ(hWJ^ = (1.®H)(1>H^) 

and 



Inll ® HH^, (30) 



C.i.o (hW,,, ^I,.,) = logdet (^I„,„ + ^^ntl ® HH^) . (31) 



3.3 Examples 

To evaluate the performance of IR and CC, several exemplary MIMO configurations are 
considered and the corresponding optimal average rate performance (p4l) fl29l) is compared 
in this section. 

We start with the case of = 1 (MISO)0, which is easy to analyze as the received signal 
expands only one spatial dimension, and hence the capacity distribution would be easy to 
obtain. For a {Lt, 1) MISO system, the capacity fl23l) becomes 

C.,.o (h^ ^I..) = log (^1 + ^ g j . (32) 

Define 

Lt 

9Lt = J2\^^\'^ (33) 

1=1 

the distribution of Cmimo (h^, ^^Ilj^ is then completely characterized by the random vari- 
able QLf If the attention is restricted to a Gaussian channel with hi ~ CA/'(0, 1), then ql^ 
follows a Chi-square distribution with 2Lt degrees of freedom and it has a PDF 

fig) = f^^^*"'e"'' 9 > 0, (34) 

where F (Lt) = t^'^^^e^^dt is the Gamma function. The CDF of git, Fcig) = Pr {git < g} - 
Io° f{'t)dt, can also be calculated and a closed-form expression is given in [28]. The CDF of 
the MISO capacity then reduces to 

.T SNR, \ ^ „1 ^„ f_ ^ (2«-l)L,l _ ^ /■(2«-l)L, 



Pr { ( , —I,, j < fl| = Pr < ^^^1 = Fa ^^^^ ) . (35) 

The equivalent IR channel capacity is given in fl32|) . and the average rate with isotropic 
Gaussian input distribution for IR can now be re-written as 

fl. = f: (fl<»> - (l - Fa (^^3^)) . (36) 



*The case of SIMO {Lt — 1) can be similarly analyzed. 
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Figure 2: The optimal average rate performance of IR and CC in a (Lj = 2, = 1) MISO 
Gaussian channel with different A^. 



As for CC, the interesting observation is that the impact of each retransmission on the overall 
capacity is only to increase the receive SNR to SNR„ = nSNR. Thus, 

C^^?UHScc)=^log(l + ^^..), (37) 



and the average rate becomes 



- - r « r._nc:MR ^c:mr • (38) 



n \ V (n- 1)SNR/ V nSNR 

n=l \ \ \ t / \ 



Unfortunately, such analytically tractable expressions are unavailable for the general MIMO 
channel, and we have resorted to numerical methods to solve the optimization problems. 
Fig. Eland Fig. [3] give the numerical resuhs of Problems (El) ([29]) for {U = 2,Lr = 1) MISO 
and {Lt = 2,Lr = 2) MIMO systems, respectively. The ergodic capacity and the optimal 
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Figure 3: The optimal average rate performance of IR and CC in a (Lj = 2, = 2) MIMO 
Gaussian channel. 



average rate of the no-feedback scheme are natural upper and lower bounds of the ARQ 
schemes [6], respectively, and they are plotted for reference. As one can see from the figures, 
IR performs better than the CC and no-feedback scheme and is close to the ergodic capacity. 
Moreover, the optimal average rate continues to increase as more ARQ rounds are allowed. 
In fact, the average rate of IR can be proved to asymptotically achieve the ergodic capacity, 
regardless of the fading distribution [6, Lemma 2]. On the other hand, the performance of 
CC is very limited, and increasing the maximum allowable ARQ rounds brings almost 
negligible gain to the optimal average rate. 
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4 Optimal LDC design: a mutual information analysis 



The previous section discusses the traditional IR and CC protocols and evaluates their per- 
formance. The IR protocol has excellent average rate thanks to two properties that allow 
full utilization of the ARQ retransmission. The first is the flexibility in rate assignment 
{Rn}n=iy which can be accomplished by either the rate-compatible channel code via punc- 
turing [4,29-34] or the recently developed rateless code [35-37]. On the other hand, the CC 
protocol has only one rate parameter R to optimize, and each retransmission can reduce the 
rate by only an integer fraction. However, one should note that the advantage of IR rate 
assignment comes at the cost of complicated packetization: each transmission could have 
very different packet length, which is undesirable in practice. From this perspective, CC is 
preferred as each retransmission is a simple repetition of the previous packet. The second 
property is that each transmission in IR is capacity-achieving, provided that the channel 
coding and space-time transmission are carefully designed. For CC, however, it is clear from 
(!26|) that such capacity optimality is not valid in general, due to the repetition nature of this 
protocol. This also contributes to its poor performance. 

This section introduces the MIMO-HARQ design based on STC, which combines the advan- 
tages of both IR and CC. The STC-based protocol has the IR capacity-achieving property 
for each retransmission, while still enjoys the repetition nature of CC. Information-theoretic 
tools will be used in this section for the STC design of MIMO-HARQ. In Section [5], this 
design problem will be revisited from the decoding error probability perspective. 

The idea can be best understood from the following simple example. Consider Lt = 2, Lr = 1 
and N = 2. For such configuration, the celebrated Alamo uti code gives 

= [xi,X2]. (39) 

The traditional use of Alamouti code would transmit X in two channel uses. With ARQ 
feedback, however, one can separately transmit the first and second column. To be more 
specific, xi is sent first. The receiver gets yi = hiSi + h2S2 + zi and jointly decodes (si, S2) in 
a maximum-likelihood (ML) manner. In case the first decoding attempt fails, a NACK will 
be sent back to the transmitter asking for the transmission of the second column X2. With 
both columns transmitted, the receiver can perform the usual Alamouti decoding to recover 
Si and S2 in the second decoding attempt. 

The idea behind this scheme is that if the channel matrix is "nice" , it may be unnecessary 
to send a full-diversity ST codeworclf]. The first transmission, which only sends xi, is an 
aggressive attempt to exploit the channel. If the random channel is not good enough, we step 
back to the Alamouti scheme, with the help of ARQ feedback. It should be mentioned that 
HARQ schemes jointly considering STC and packet retransmission has been studied in the 
literature, e.g., [15,17], but none of the designs are from an information-theoretic perspective 

''This is further discussed at the end of Sec. O 
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and hence do not give the best performance one can hope from STC-based HARQ. On the 
other hand, the information-theoretic studies of such schemes [8, 10] are mostly focused on 
the high SNR diversity-multiplexing-delay tradeoff and do not study the optimal STC design 
for a finite SNR. 



4.1 General framework 



The Linear Dispersion Code (LDC), proposed by Hassibi and Hochwald [38], is a general 
STC design that incorporates most of the known STBCs as special cases. This tool is used 
and a LDC framework is formulated to study the performance of STCs for the MIMO-HARQ 
design. Consider a {Lt,Lr) MIMO system. For each LDC codeword K modulated complex 
symbols are transmitted, and T is used to denote the maximum number of time intervals 
of the LDC codeword for all ARQ rounds, i.e., T = '^n=i'^n- A LDC codeword can be 
expressed ClS db linear combination of the symbol 



K 



k=l 
K 

= ^SfcCfc + 4Dfc, (40) 

k=l 

where Sk = Oik + jPk are the modulated symbols, and A^, B^, C^, G C^*^"^ are the LDC 
spreading matrices with 

Afc = Cfc + Dfc, 

Bk = Cfc-Dfc. (41) 



For the HARQ transmission with a maximum N rounds, the overall LDC codeword X shall 
be divided into N sub-codewords 

X= [Xi,X2,--- ,X^], (42) 

where X„ G C^*^^" is for transmission at ARQ round n. The LDC spreading matrices can 
be similarly divided 

Afc = [Afc^i, Afc^2, ■ ■ ■ ; Afc^Tv] , 

Bfc = [Bfc^i, Bfc^2, ■ ■ ■ , Bfc^Ar] , (43) 

or 

Ck = [Cfe_i, Cfc,2, ■ ■ ■ 5 Cfc^Ar] , 

C,. = [D,,i,Dfc,2,--- ,Dfc,^], (44) 

^As our focus is to study the mutual information of the LDC structure with HARQ, the effect of channel 
coding length L„ is irrelevant and is ignored in the discussion. 
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where A^^^, Bfe^„, Ck,n, Dfc,„ G C^*=^'^". Defining 



rp{n) 








= [Afc_i, Afc^2, ■ 


; ^k,n 




= [Bfc,i, Bfc^2, ■ ■ 


■ ■ , Bfc^ri] 




= [Cfc,i, Ca;,2, ■ ' 






= [D/fc^i, Dfc^25 ■ 


5 D/j^fl 



(45) 



(46) 



the accumulated receive signal after the n-th transmission can be written as 



Y(") = 7^HX(") + Z("), (47) 



with 

x(") = [Xi,X2,---,X 



X:a,A(")+j/3,Bi") (48) 

k=l 

f^s.CP + slBtl (49) 

k=l 



The capacity of the signal model P7|) is denoted as C^^\li), which is determined by the LD 
structure (HSjl or fjl^ . Meanwhile, notice that the effective rate for ARQ round n is 

^^"^ = (50) 

where R is the channel coding rate, which is a design parameter. Now directly applying the 
general result (fT6ll . the average rate of LDC-based MIMO-HARQ can be expressed as 

N 

= ^(i?(")-i?("+i))Pr|4"^(H)>i?(")} 

ra=l 

n=l *~ ^ 

where we use the convention T^^~^^^ = oo such that = 0, and T/v+i/T'^^"'"^) = 1. 

The general problem of average rate maximization for LDC-based MIMO-HARQ can be 
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formally casted as 



maximize 
subject to 



JV 

n = l 



i? > 



Tn is a positive integer 



(52) 



Tr(Ef=i (A,,„Ag„ + B,,„Bj„)) =2L,T„,Vn = l,--- ,iV. 



Notice that C^^ (H) implicitly relies on {Aj, Bj, Tj}"^-^. The power constraint in Problem 
(132]) will be discussed in more detail in Section I4.4[ 



4.2 Design criterion 

Directly solving Problem (|52i) . either analytically or numerically, is difficult, but there are 
some key properties that will help guide the optimal LDC design. The first observation is 
that the average rate is monotonic with the LDC capacity, and an achievable upper bound 
is the capacity of the MIMO channel. This leads to the following design criterion for the 
LDC design of MIMO-HARQ. 

Criterion 1 (Capacity-based LDC design criterion) The LDC should satisfy 



for any ARQ round n = 1, ■ ■ ■ ,N. 

Notice that the intuition behind this criterion is that each LDC transmission should create 
an equivalent MIMO channel that is capacity lossless. Under the constraint of isotropic 
Gaussian distribution, this criterion gives the largest capacity for each transmission, which 
then minimizes the probability of decoding failure for any given R and {T„}^^^. 

With Criterion [H the original Problem (1521) is greatly simplified, as the dependence on 
{A„, Bn}n=i disappears. The problem becomes how to find the optimal R and {T„}^^-^ 
that maximizes the average rate. In most of the remaining work a special case of T„ = 1, 
Vn = 1, ■ ■ ■ , is considered. There are two reasons to focus on this special situation. 

1) Optimizing over all possible combinations of is complex even for moderate N. 

This requires solving a single- variable optimization problem for each possible combina- 
tion of {Tn}^^i satisfying J2n=i'^n = T and T„ being a positive integer. The optimal 
solution may change with respect to the operating SNR or the channel distribution. 
In practice {Tn}^^^ may be pre-determined with the choice of LDC, which leaves only 
one variable R to optimize over. 




(53) 
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2) Choosing T„ = 1 ensures that each ARQ round is of minimum delay, which is an 
important advantage. Meanwhile, this minimum delay may also result in the best 
throughput, as it is empirically observed as a good balance between ^(i^^l^+i) and 

Pr jCiJ"^ (H) > Intuitively, large T„ will result in a marginal gain in the decoding 

successful probability, but a significant decrease of the multiplicative coefficient. 

To summarize, we give the following corollary. 



Corollary 1 Assume isotropic Gaussian input distribution over all available transmit an- 
tennas. For a given set o/{T„}^^-^, the optimal average rate of a LDC-based MIMO-HARQ 
protocol is 



N 



R\d = max 



RTn+l 



SNR, 



R 



In the special case ofT^ = 1, Vn = 1, ■ ■ ■ , A^, the optimal average rate becomes 
^ Prlc^..^Ju.^l.]> 

n N 



(54) 



(55) 



Criterion [T] and Corollary [T] are based on the mutual information of the equivalent MIMO 
channel. To more conveniently evaluate the optimality of any given LDC structure, it is help- 
ful to develop a criterion that explicitly relies on the LDC spreading matrices {0^,0^}^^^, 
which lead to Theorem [TJ 



Theorem 1 Assume Lr > Lt and K = LtT. Define 



vec C\ 



,vec{ CP 



vec I Dl 



vec i D 



v{n) 

'k 



and 



(56) 
(57) 



jj{n) -yin) 

conj (V(")) conj (U(")) 

The LDC-based HARQ protocol {Ck,n, ^k,n}n=ik=i leads to the optimal average rate perfor- 
mance if and only if 

F(")F(")'' = (58) 

for all n = 1, ■ ■ ■ ,N. 



The proof follows almost directly from [39, Theorem 5.3.1]. The challenge is that Theorem 
[U requires to consider K > LtT^"^^ for each n = 1, ■ ■ ■ , — 1, while [39, Theorem 5.3.1] 
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only holds for K = LtT. However, a careful study of the proof of [39, Theorem 5.3.1] shows 
that F(")F(") = ^2^(n) would be necessary and sufficient for the LDC to be average rate 
optimal. 

Theorem [1] considers the general LDCs that allow the conjugation operation. This is a key 
feature in some codes, e.g.. Orthogonal STBC, but it has been argued [40] that there is not 
always a significant gain over complex LDCs, in which no conjugation operation is allowed. 
Hence, it is instructive to give the design criterion for complex LDCs, which follows directly 
from Theorem [H 

Corollary 2 For a complex LDC 

x(-)=x:..cr) (59) 

k=l 

with Lr > Lt and K = LtT, it is average rate optimal if and only if 

U(n)uW^ = I^^^(„,, (60) 

for all n = 1, ■ ■ ■ , N . 



4.3 Evaluating existing LDCs 

Previous sections analyzed the optimal average rate performance of LDC-based MIMO- 
HARQ and proposed the design criterion of the corresponding space-time structure. This 
section will attempt to use this general framework to evaluate existing LDC structures and 
their usage in the MIMO-HARQ setting. Although Theorem [1] and Corollary [2] provide 
straightforward means to evaluate given LDCs, we shall proceed to study the equivalent 
MIMO channel capacity and use Criterion [1] whenever it is feasible and simple. Directly 
evaluating the equivalent MIMO channel capacity and comparing to the physical channel 
capacity shed light on the (sub)optimality of the LDC. We will start with the simple MISO 
Lf = 2, Lr = 1, N = 2 setting, in which analytical (sub)optimality can be rigorously shown. 
For the more complicated MIMO settings, numerical simulations are performed as a main 
tool to evaluate existing LDCs. 

4.3.1 MISO Lt = 2,Lr = l,N = 2 

This is the simplest MISO setting, in which the celebrated Alamouti code was invented [41]. 
The following LDC-based HARQ protocols are studied. 

Spatial Multiplexing with Repetition: 
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This is the multiple-antenna version of the CC scheme. Spatial multiplexing vector x = 
[x(l),x(2)]* is repeated upon a NACK, and the overall codeword after N = 2 transmissions 



IS 



x(l) x(l) 
x(2) x(2) 



Since ARQ round 1 is a spatial multiplexing, it is capacity-optimal: 

SNR 



(61) 



C«(h*) = C^i.o([h(l) h(2)] 



I 



Lt 



, , , SNR 

log 1 + -T^92 



(62) 



where g2 is defined in fl5^ . The simple repetition in ARQ round 2, however, leads to a 
capacity loss: 



h(l) h(2) 
h(l) h(2) 



SNR 



Lt 



-log(l + SNR(72) 



-C^i^o (h*,SNRIiJ 



< a 



(63) 



for SNR > 0, which proves the strict suboptimality of this protocol. 
Antenna Switching (AS): 

This scheme activates only one antenna at a time, and rotates the active antenna to achieve 
spatial diversity. The overall codeword is 



X 
X 



(64) 



This is apparently a low-rate suboptimal scheme with 



log(l + SNR(?i)<a 

1 



iog(i + sNR(72) <a 



and hence neither transmission is capacity-optimal. 



(65) 
(66) 



Alamouti Code: 



19 



The overall codeword after N = 2 transmissions is 



x(l) -x(2)' 
x(2) x(l)* 



(67) 



The first transmission is a spatial multiplexing, which is optimal. It is also well known that 
the overall Alamouti structure also achieves the channel capacity [38]. Hence, 



C 



(1) 



Alamouti 



(h*) — C^ig^outi (h*) — Cmimo 



t SNR 
- — 1/ 



(68) 



and this protocol satisfies Criterion [TJ 
Cyclic Delay Diversity (CDD): 

CDD was proposed in an OFDM setting [42] but the code itself presents an interesting 
space-time structure. It cyclically rotates the previously transmit column such that each 
symbol is sent from different antennas in different times slots. The overall codeword in the 
{Lt = 2,Lr = 1) MISO is 

"x(l) x(2) 
x(2) x(l) 

ARQ round 1 again is a spatial multiplexing, and hence is optimal. For ARQ round 2, the 
equivalent channel matrix is 



(69) 



H 



eq, cdd 



h(l) h(2) 
h(2) h(l) 



(70) 



and 



c'2 (h') 



^logdet (I2 + ^^Heq,cddH2,cdd 

1 + ^92 SNR Re{h(l)h(2)*} 

SNR Re{h(l)h(2)*} . , 2 

1. SNR, ^.,2\ 1. SNR 



- log det 
2 ^ 



i + -^92 



log 1 + 
2^1 2 



SNR 



< log(^l + 

— C ■ I h 



2 ^ 
t SNR 



|h(l) + h(2)|^J +-log( 1 + 
|h(l)p + |h(2)p) 



hfi)-hm 



(71) 



This suggests that the CDD-based protocol is suboptimal. 
Numerical Simulations: 

In addition to the previous analysis, numerical optimization is performed to obtain the 
average rates of several LDCs and the results are reported in Fig. HI Optimal LDC, which 
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Figure 4: The optimal average rate performance of several existing LDCs in a (Lj = 2,Lr = 
l,N = 2) Gaussian MISO-HARQ channel. 



satisfies Criterion [T], and IR are plotted as references. It is interesting to observe that 
although both AS and CDD are analytically proved to be sub-optimal, the average rate 
performance of CDD is very close to the optimum (almost negligible performance loss), 
while AS is extremely sub-optimal. This is no surprise as CDD utilizes the spatial degrees 
of freedom in a more efficient way than AS. SM-based repetition gives the performance of 
CC and is not far away from the optimal LDC-based protocol. Alamouti, as predicted in 
the analysis, is the optimal LDC. 



4.3.2 Full MIMO 



Analytically proving (sub)optimality becomes infeasible in most of the full MIMO settings, 
as the comparison of mutual information of a given LDC {A„, B„}„^i to the MIMO channel 
capacity is difficult to make even for moderate N. This section thus focuses on numerical 
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Figure 5: The optimal average rate performance of several existing LDCs in a (Lj = 2, = 
2,N = 2) Gaussian MIMO-ARQ channel. 



simulations to get insight on existing LDCs. 



Fig. [5] plots three well-known LDCs in a {U = 2, = 2) MIMO channel with N = 2: the 
Golden Code (Golden) [43], Damen, Tewfik and Belfiore's code (DTB) [44], and Hassibi and 
Hochwald's LD code (HH) [38, Equation (31)]. These three LDCs are known in the literature 
to be capacity optimal, and in fact it can be proveiillthat they satisfy Corollary [2] or Theorem 
[H and is average rate optimal as well (the first column of the codeword is capacity lossless). 
Similar to the mutual information maximization approach of LDC design [38], typically there 
exist many LDCs that satisfy Criterion 1. In practice, one can search among these codes to 
further consider other criteria [40], such as the diversity and coding gain. 

For a {Lt = 4, = 2) MIMO-HARQ system, Figs.[6]and[7]give the simulation results for = 
2 and = 4, respectively. The LDCs under consideration are: Orthogonal STBC (OSTBC) 

^This can be done by direcly verifying Corollary [2] or Theorem [1] using their corresponding LDC spreading 
matrices. 
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10 15 20 25 30 35 



SNR [dB] 

Figure 6: The optimal average rate performance of several existing LDCs in a (Lj = 4, L,, = 
2, = 2) Gaussian MIMO-ARQ channel. 

[45, Chapter 7.6], Diveristy Embedded STC (DE) [46], ABBA code [47], Double ABBA 
code (DABBA) [48], Jafarkhani's quasi-orthogonal STBC (J-QOSTBC) [49], Papadias and 
Foschini's QOSTBC (PF-QOSTBC) [50], Ran, Hou and Lee's QOSTBC (RHL-QOSTBC) 
[51], and Double Space-Time Transmit Diversity (DSTTD) [52]. CC, IR and optimal LDC 
are again plotted as reference curves. 

Most existing LDCs for Lf = 4 are obtained using Lf = 2 LDCs as building blocks, which 
use 2 time slots. Hence, unlike other examples in this paper, the HARQ protocols for N = 2 
use Ti = T2 = 2 and the results are reported in Fig. O The first observation is that even 
the optimal LDC has an average rate that is very close to CC. This suggests that with 
this specific configuration, one should not expect much average rate gain from LDC over 
simple repetition protocol. However, a well designed LDC-based protocol may have other 
advantages, such as diversity gain and simple decoding. The second conclusion is that among 
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Figure 7: The optimal average rate performance of several existing LDCs in a (Lj = 4, = 
2, iV = 4) Gaussian MIMO-ARQ channel. 
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all the LDCs, only DSTTD performs close to the optimal LDQj. Among the other LDCs, 
five (ABBA, DABBA, DE, J-QOSTBC and PF-QOSTBC) have very similar optimal average 
rate, while RHL-QOSTBC is worse and OSTBC gives the worst performance. The extreme 
suboptimality of OSTBC is again due to its low-rate property (rate 3/4). 

The conclusions are different in the case of = 4, in which each ARQ round sends one 
column of the overall ST codeword matrix. From Fig. [TJ there is a notable gain of optimal 
LDC over CC. Meanwhile, none of the above LDC^ are optimal, which calls for the search 
of novel LDCs that approach the performance of optimal LDC. 



4.4 Remarks 

Remark 1 Power Constraint 



We would like to constrain the long-term average transmit power in Problem (1521) . Two 
different types of power constraints can be considered. 

1) Each ARQ round has the same power. 

In this constraint, regardless of how many transmissions actually take place, the overall 
transmit power is constant due to the equal power allocation over ARQ rounds. The 
power constraint in the original problem (152!) falls into this category: 

Tr ||fj(A,,„A^„ + B,,Xjj =2L,r„,Vn = l,---,Ar. (72) 

Note that with this constraint, each transmit antenna or each time slot can have dif- 
ferent transmit power. Also the real and imaginary parts of the transmit constellation 
can have unequal power. 

Similar to [38], more stringent power constraint can be posed to replace (1721) : 

Tr (A,,„A^„) = Tr (B,,.BJ„) = (73) 

or 

T 

Afc.nAfc „ = Bfc,„B;, „ = -^Il„ (74) 

where ( I73p ensures real and imaginary parts are transmitted with the same power, while 
( in|) further forces energy to be spread equally in all spatial and temporal dimensions. 
These power constraints can be incorporated into Problem (l52l) . 



^It should be noted that DSTTD is in fact capacity suboptimal. 
^DSTTD only uses two time slots and thus cannot be used with = 4. 
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2) Dynamic power allocation among ARQ rounds. 

Instead of always having a constant power for each ARQ round, one can further allow 
power allocation among the N ARQ rounds. For example, giving more power to ARQ 
round 1 will increase the probability that the first transmission succeeds. However 
this comes at the price of decreasing the power of potential ARQ rounds 2 to N, if 
the overall average transmit power is kept constant. There is a serious challenge to 
maintain a constant average power. The reason is that except for the first transmission, 
ARQ rounds 2 to A happen only with a non-one probability. This challenge has been 
addressed in the scalar channel case [6], and the solution is to average the power 
consumption with respect to the channel fading distribution. 

Assuming that ARQ round n has a total transmit power p„, the actually consumed 
power p is a random variable with PMF 

J27=i pi'^i if A ~A A . vjrn — 1 AT — ^ 
^^^Vt — (75) 

Note that dynamic power allocation also affects the equivalent MIMO channel capacity 
after ARQ round n: C|J^(H) is a TDMA-type one where different portions of the code 
have different power [6]. 

To maintain an average total transmit power SNR, (pi, ■ ■ ■ ,Pn) should satisfy 

E %^^P^ { A , ■ ■ ■ , AIT, A } + ^^K^Pr { A, ■ ■ ■ , = 1. (76) 



Remark 2 Chase, LDC, and IR 



Let us revisit the average rate for IR (|24l) . CC ( |29l) . and LDC (l52ll . It is straightforward to 
verify that if the rate assignment in IR is 

R 

= t(") is a positive integer, Vn = 1, ■ ■ ■ , A, (77) 

the average rate expression of IR and LDC will be the same, provided that Criterion [T] is 
satisfied. This demonstrates that the advantages of IR over LDC is the flexibility in rate 
assignment for each ARQ round, which comes at the cost of unequal packetization. 

In reality, nevertheless, the freedom in rate assignment of IR is very limited due to practical 
limitations. Considerations such as packetization or slotted multi-user transmission usually 
ask for constant-length retransmissions, which leads to the same IR rate assignment as in 
LDC. With this constraint, the optimal average rates of IR and LDC are the same. In such 
cases, LDC has the potential advantage of receiver complexity over IR. Recall that IR requires 
code combining at the receiver - punctured codeword symbols need to be combined with the 



26 



previous transmissions for another decoding. LDC-based protocol, on the other hand, only 
needs packet-level combining. Thus, the receiver complexity is shifted from channel decoding 
to front-end demodulation/detection. If the LDC is well designed, such detection complexity 
can be moderate, e.g., the Alamouti code. This can be an important advantage of LDC. 

At the same time, LDC can also be viewed as an enhanced CC scheme. The reason is that 
in the LDC-based protocol each retransmission is essentially just a repetition of previously 
transmitted symbols - no new information is sent. However, the "repetition" in LDC is 
handled in a smart way such that the equivalent channel capacity is increased to approach 
the MIMO channel capacity, while the simple repetition CC scheme suffers from capacity loss. 
This can also be explained from the diversity perspective, which will be further developed 
in Sec. [51 The CC protocol naively repeats the same codeword upon each NACK, which 
cannot efficiently exploit the spatial diversity. A good LDC-based HARQ protocol spreads 
out the symbols onto different antennas upon each retransmission request to better exploit 
the spatial diversity. This becomes especially important in a quasi-static fading channel, as 
there is no time diversity to exploit. 

Remark 3 Switching Between Multiplexing and Diversity 

Compared to the STC-based HARQ scheme, there is a similar work [53] on switching be- 
tween multiplexing and diversity, which is worth some comments. It is easy to see that 
the Alamouti example essentially performs a switching between multiplexing and diversity, 
just as [53] does. However, one can see that the ARQ-feedback-based approach is supe- 
rior to the method in [53]. With ARQ feedback one actually sets the default choice to be 
full multiplexing. Hence, data communication takes place without even knowing whether 
one should use multiplexing or diversity. Switching to the diversity scheme only happens 
whenever necessary. This is possible because full multiplexing is embedded in any diversity 
schemes, which is not exploited in [53]. For the scheme in [53], the transmitter needs to 
wait for the feedback to inform which scheme to use before the data communication can take 
place. Another advantage of our approach is that ARQ feedback informs the decoding status 
to the transmitter, and the decision of switching to diversity scheme is made only when the 
receiver fails to decode the full multiplexing transmission. This is better than the selection 
criteria of [53], which use the minimum Euclidean distances and Demmel condition number 
of the MIMO channel to make the decision. Notice that the ultimate goal is to improve the 
decoding error performance, and hence decoding failure should be the best decision metric 
to determine switching. Our approach directly relies on the decoding error event, while [53] 
uses some indirect performance measures, which can only approximate the decoding error 
event. 
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5 Optimal LDC design: an error probability analysis 



Space-time codes design has followed two related but different paths. The original approach 
is from the viewpoint of detection theory. The focus is on studying the pariwise error 
probability (PWEP) [54-57], which has led to the well-known rank and determinant criteria 
that use diversity and coding gain to compare different STCs. A different approach is from 
the standpoint of information theory. It views STC as a modulation and study the optimal 
structure that preserves the channel mutual information [38,58]. There are also works that 
combine these two views [39,53]. The previous section investigated the optimal LDC design 
from an information-theoretic point of view. In this section, an error probability analysis of 
the optimal LDC-based HARQ protocol is performed, and the corresponding design criterion 
is presented. 



5.1 n-th pairwise error probability 



Assume that the signal vector s = [si, ■ ■ ■ , sxf is chosen from a uniformly distributed set S 
with cardinality M: 

5 = {so,--- ,sm-i}. (78) 
The probability of a decoding error after the n-th ARQ round can be written as 



A/-1 



Pi") = {^^'^ ^ - - - , s(") ^ s, |s, was sent} 



i=o 



M 



AI-l 



M-1 M-1 



U---U {c<0'---'Qa<o}iH,s, 



j=0 

M-1 M-1 M-1 



n=0 in=0 



(79) 



^0) 



'1) 



j=0 ii=0 in=0 



where ( IHTj) comes from the traditional union bound [59], Sj is the transmitted vector, s^"-* 
is the detected vector after ARQ round n, Qf^j is the pairwise decision metric in the n-th 
ARQ round, which is defined in the following. The optimum ML decoding rule for deciding 
between two possible codewords X*^"')(sj) and X'^"')(sj) for a given channel realization H and 
the receiver observation Y^") is 



Define the ML metric 



^HX(")(s.)||| ; IIyW 



SNR 



Lt 



HXW(s,)||^ 



SNR 



HX(")( 



(82) 



(83) 
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Figure 8: An example illustrating different error events in different ARQ rounds. Sj (red) is 
transmitted, and s*^") is the ML detector output at ARQ round n, which are all erroneous 
(blue). 

and qI^j is defined as 

Qin) ^y(n) _yin)^ (84) 

For n = 1, (IHTl) reduces to the traditional union bound for STC: 

Af-l M-1 

^ M E E ^ < 0|H, s,}] . (85) 

j=0 i=0 



For Rayleigh fading, 



Pr{g,,<0|H,s,} = g| (86) 
where d'^^{i,j) is the squared Euclidean distance between two received codewords Y^^"'' and 



y\'"\ with Yf ^ = ^^HX(")(si). Tools such as the Chernoff bound [55,56] or the tighter 

bound in [57] can be applied. Meanwhile, assuming high SNR would further simplify the 
upper bound and it has led to the rank and determinant criteria. 

The performance analysis of MIMO-HARQ, however, is much more involved. The first 
complication comes from ( jSTl) . Unlike STC where only the PWEP is needed, analyzing 
MIMO-HARQ requires the n-th pairwise error probability (n-PWEP): 

Pr{Q«.<0,---,Q5<0|H,s,}, (87) 

in which the decision metrics j | are correlated random variables. The complication 

is also intuitively understandable. The error event under consideration is when all the n 
decoding attempts are erroneous. However, the decoding errors in different ARQ rounds 
could be very different (especially for short codes), depending on the channel realization and 
instantaneous noise. Fig. [H] is one such example. 
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To analyze the n-PWEP, let us start with the statistics of Q^^^^ . To simphfy the derivation, 
define 



D 



(n) 



SNR 



H(X(")(s,)-XW(s,)), 



(88) 



and 



(n) 


= vec 


(y<">) , 




= vec 






= vec 


(Z<"') . 



(n) (n) 
y) : 



(89) 



Now conditioning on H and assuming is transmitted, we have 

Qg) = ||D('!) + Z(")||i-||Z(")||2. 



IDS)||^ + 2Re{lV (Dg.)zW^)} 



(n) 



(90) 
(91) 



where (i, j)^ = | 1 1|, is the squared Euclidean distance, and = 2Re |Tr [pf) Z^"^ 
It is easy to see that w\^} is a real Gaussian random variable 



Wl""} = TrfDi".)z(")'' + Di".)''z(")) 



with 



0, 
E 

E 



(92) 



(93) 



(94) 



which uses the fact that Z^^^ has complex circularly symmetric Gaussian entries with unit 
variance. 



The n-PWEP now becomes 

Pr {q5 < 0, • • • , q5 < 0|H,s,} = Pr {<] < 
Notice that 



,w,s<-<^(^n,jriH,s. 

(95) 
(96) 



w 



(n) 



(1) 



(n) 
in,j 
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is an n-dimensional real Gaussian random vector. Thus, the only remaining problem is to 
obtain the statistics of w'^"^. It is easy to get the mean 



E[w(")] = 0. 

As for the covariance matrix Efw^")^], the focus is on 

R^(„,(fc,/) = E[iy(^>5],V/>fc 
due to its symmetric property. The important observation is that 



(97) 



(98) 



.(0 



vec (Z«) 

vec([zW,Z;_fe]) 



(99) 



and 



i(0 

it,j 



.(0 



(k) (k) 

H,3 



Hi,j,l-k 



(100) 



Hence 



R^(„)(A;,0 = E 



"d;:>+d«/zw)(z(')''da+d£/z<')) 



2Re{E[d<:-V>(d«;;/z«)']}, 



(101) 



which again uses the circularly symmetric property of the i.i.d complex Gaussian random 
matrix Z. The following derivation is straightforward: 



E 



dSV)(d<;;/z<'))' 



E 
E 



dW/z<')(d£''z« + d..„,_,%_.)' 

d«/z«(d;3''z<'))' 



(102) 



H 



2(fc)z(fc) 



H 



where for matrices A and B of the same dimension, 

(A, B)^ = Tr (AB^) 



(103) 



(104) 
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is the Frohenius inner product [60], which induces the Frobenius norm. Equation (11021) is 
because z^'^) and z;_fc are independent. 



Now the covariance matrix Rw(n) is fully characterized as 

R^(„,(A;,/) = 2Re{(D5,D5)^},V/>fc 

R-wM = R^(n) {l,k),'^l < k, 



and we conclude 



w(") ~ AfiO, R, 
The n-dimensional joint Gaussian PDF of w*-"-* is 

1 



,{n) , 



/w(") (w) 



(27r)"/2^/detR; 



: exp 



(105) 

(106) 
(107) 



:i08) 



(109) 



In the existing literature of Gaussian tail distribution, much attention has been paid to 
either one- or two-dimensional Q functions [28]. It is clear that to analyze MIMO-HARQ 
with deadline A^, one needs a general definition of the n- dimensional Q function for a real 
Gaussian random vector x ~ A/'(0, Rx) (assuming zero mean without loss of generality) 



Qn (x, Rx) = 



(27r)'^/2ydetR; 



poo 






■ / exp ( 


/x(l) 


Jx(n) \ 



;iio) 



Finally, the union bound on Pi"^ in (ISlll becomes 

^ M-l M-l M-l 



j=0 ii=0 in=0 



where 



(111) 



(112) 



5.2 Design criterion 



The previous section studies the detection error probability after the n-th ARQ round. Using 
the traditional union bound shifts the focus to the n-PWEP. With the help of ra-dimensional 
Q function and the distribution of the decision metrics, the n-PWEP is successfully cal- 
culated. However, the resulting union bound expression is complex and does not provide 
insight into the optimal LDC design. This section is devoted to developing the LDC-based 
HARQ design criterion from the detection error probability point of view. More specifically. 
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the focus in on studying the diversity order of LDC-based HARQ0. The main result is stated 
in the following theorem. 

Theorem 2 The optimum diversity order for any LDC-based HARQ protocol satisfies 

div; = L,min{T("),Lt} (113) 
for any ARQ round n, \fn = 1, ■ ■ ■ , N . 

Proof: The following upper bound is straightforward: 

Pr{A,--- ,A} <Pr{X}. (114) 
Using this bound, the error probability ( 1791) can be bounded as 

^ M-l 

Pi") < — ^ Pr {s(") ^ Sj\sj was sent} . (115) 
i=o 

The right hand side of flllSp is in fact the detection error probability of a STC with codeword 
X = X*^") G C^*^^^"\ The achievable diversity order for such STC is well-known to be 
L^min{T("),Lj. Hence 

div,; > L^min{T("),Lj. (116) 

On the other hand, the fundamental diversity order offered by the MIMO-ARQ channel has 
been established in [8,9], which leads to 

div; < L,min{T("),Lj. (117) 

Combining (11161) and (I117P completes the proof. 



Theorem [2] leads to the following diversity-driven LDC design criterion for MIMO-HARQ. 



Criterion 2 (Error-probability-based LDC design criterion) The LDC should satisfy 
that after the n-th ARQ round, = 1, ■ ■ ■ , N , the optimum detector has a diversity order 
Lr min{T("),Lj. 



^°In this work, the diversity order is defined as 

. logPe(SNR) 
div = — hm — - — — — — 
SNR^oo logSNR 

for a fixed-rate STC. This is the traditional definition of diversity, which is different from the one defined 
in [7] and many following DMT papers. 



33 



This design criterion has some practical advantages over the capacity-based Criterion [H as 
it is quite intuitive and easy to check. Meanwhile, it is also observed that those LDC-based 
HARQ protocols that satisfy Criterion [T] also tend to satisfy this diversity criterion. 

Traditional STC design focuses on both the diversity and coding gain. A typical design 
would use full diversity as a constraint and then search for good coding gain. Our effort has 
been only on diversity due to the following two reasons. The first is that coding gain for 
MIMO-HARQ is not easy to define or compute, due to the complication in the n-PWEP. 
Theoretically, each ARQ round is equivalent to a STC X = X*^"^ G C^^^^''"\ for which the 
coding gain can be computed. It is unclear how to optimize these N correlated coding gains 
simultaneously. The second reason is that coding gain can be obtained from outer channel 
code as well. The diversity gain, however, has to be exploited by the proper LDC design. 
This is because the quasi-static channel is considered, and the only source of diversity is 
from the spatial domain. 



5.3 Numerical Examples and Discussions 

To gain more insight into the LDC-HARQ design based on the error probability analysis, we 
return to the Alamo uti example ( l39l) for N = 2. Studying this simple example leads some 
interesting observations, which reveals the significance of ARQ feedback in different settings. 

A software simulation system for {Lt = 2, = 1) MISO i.i.d. Gaussian quasi-static fad- 
ing channel with maximum ARQ rounds = 2 is built. Alamouti-based protocol (!39l) is 
implemented, as well as the traditional repetition scheme (CC). Maximum-likelihood de- 
coding is performed for both ARQ rounds, which gives the optimal decoding performance. 
We consider both uncoded and coded transmission with a Gray-coded QPSK constellation 
on each antenna. For the coded case, the standard 64-state rate-1/2 binary convolutional 
codj"1 with octal generators (133, 171) is implemented at the transmitter, which is followed 
by a standard row-in-column-out block interleaver. Each coded packet consists of 100 sym- 
bols. The receiver implements ML soft-output detectors with exact bit log-likelihood ratio 
(LLR) computation concatenated with a bit deinterleaver and a soft-input Viterbi algorithm 
convolutional decoder. Finally, it should be emphasized that as discussed in [6], the ran- 
dom channel generation in the simulation guarantees that the empirical channel distribution 
matches the true one. 

Fig.[n]reports both the packet error rate and the average rate performance as a function of the 
average receive SNR. As has been analytically predicted in this work, the traditional ARQ 
approach which simply repeats the previous packet is significantly inferior to the Alamouti- 
based protocol in all settings. There is a clear diversity loss for the CC protocol, which also 
results in the degraded average rate performance. It should be remarked that the gains of 
Alamouti-based HARQ over CC are not only in packet error rate and average rate, but also 

^^The rate loss due to the termination bits is ignored. 
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in the receiver soft-output ML detection complexity. Alamouti detector orthogonalizes the 
two transmitted symbols for independent bit LLR generation, while the CC protocol does 
not render such simple detector but an equivalent 2x2 MIMO system has to be formed to 
generate bit LLR, which is of higher complexity. 

In general, the gains of HARQ protocols over direct transmission without ARQ feedback are 
two fold: 

1) throughput; 

2) decoding error probability. 



The interesting observation from Fig. [H] is that these gains behave quite differently in different 
system configurations. The average rate gain for using ARQ feedback is remarkable for all 
considered protocolf^, as can be seen in Fig. 9(b) However, the packet error rate gain of 



the two considered HARQ protocols over no-ARQ transmission is distinguishable for the 
uncoded case, but is almost negligible for the coded one, even if the channel coding is not 
capacity-approaching and the block length is short. This can be explained from a joint 
consideration of the mutual information analysis in Sec. H] and the error probability analysis 
in Sec. El Let us take iV = 2 as an example. The decoding error probability of using ARQ is 



as compared to the no-ARQ case 



Hence, the gain of HARQ is 



Parq =Pr{A,A}, (US) 



Pno-arq =Pr{A}. (119) 



Pno-arq-Parq = Pr { A } " Pi { ^1 , ^2 } 

= Pr{^i,:4^} (120) 
> 0. 

Notice that the event {^1,^2} means the receiver has a correct decoding in the first ARQ 
round, but it makes an incorrect decoding as it performs decoding based on both packets. 
This event is possible in the uncoded there is a good probability that the noise 

realization in round 1 is small (hence the decoding is correct), but very large in the second 
round such that the decoding based on both transmissions fail. Recall that the channel 
transfer matrix keeps constant over ARQ rounds and hence the only source for this to 
happen is the noise realization. This explains the gain in the uncoded case. However, as the 
channel coding is implemented and the packet length is increased, the impact of finite noise 
realizations will eventually vanish. In fact, for the capacity- approaching coded transmission 



^^Some uncoded average rates seem to be better than coded ones in low SNR. This is because the coded 
packet is longer than the uncoded one, which counteracts the benefit from (short) channel coding. 
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SNR [dB] 

(a) PER 




Figure 9: Performance comparison of two HARQ protocols (Alamouti and CC, N = 2) and 
no ARQ, with both uncoded and BICM coded transmissions in a [Lf = 2, = 1) MISO 
Gaussian channel Fig. 9(a): Packet Error Rate (PER) versus SNR. Fig. 9(b) Average Rate 
(bits per channel use) versus SNR. 
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as discussed in Sec. [3] and HI the noise realizations will be totally averaged out, which leads 
to Equation (jS]) and hence 

Pr{^i,:4;} = 0. (121) 

In summary, even with a relatively strong channel code with reasonably long block length, 
the two information-theoretic observations in Sec. [3] are generally valid, and the ARQ re- 
transmission leads to only an average rate gain but very little error probability gain. On the 
other hand, very short coding length makes the ARQ retranmission even more beneficial, as 
the decoding error probability is improved prominently in addition to the average rate gain. 

One might notice that there is only one degree of freedom for a {Lt = 2, = 1) MISO 
channel, but two independent symbols are sent in the first ARQ round. This observation 
reveals the essence of the HARQ gains. The concept of degrees of freedom is from the 
ergodic capacity, which averages over the channel distribution. For a certain subset of all 
possible channel realizations (together with appropriate noise realizations if the finite noise 
impacts the decoding), sending two independent symbols per channel use can result in a 
successful transmission. In this case, the spectral efficiency is doubled compared to the no- 
ARQ transmission, which is refiected in the average rate gain. If the channel realization is 
"bad" , the second round transmission switches the entire protocol back to the normal STC 
transmission. 



6 Conclusions 



This paper studies the fundamental performance of Hybrid ARQ protocols in a multiple- 
antenna channel. Unlike previous works that focus on the high SNR asymptotics, a general 
framework to optimize the average rate for any HARQ protocol with fixed SNR is presented. 
This general result is then applied to study two well-adopted HARQ protocols, incremental 
redundancy and packet combining (Chase). Space-time coded HARQ transmission is inves- 
tigated under this framework, in which the capacity-based LDC design criterion is derived. 
Several existing LDCs are evaluated with this criterion, both analytically and numerically. 
A different design criterion based on the error probability analysis of LDC-HARQ is also 
presented. Numerical examples reveal different types of advantages of LDC-HARQ, and 
how they vary in different settings. 

There are several interesting problems that have not been addressed in this work, which 
are the subjects of potential future work. For several MIMO and ARQ configurations, the 
existing LDCs are proved to be optimal in terms of average rate. However, there are some 
other settings where none of the known codes approach the optimal performance, e.g., the 
{Lt = 4, Lj. = 2, N = 4) configuration in Fig. [71 Numerical design/search of optimal codes 
(e.g., similar to the approach in [38]) is an interesting research topic. Another possible 
direction is to incorporate the receiver decoding complexity into the LDC design criteria. 
Recall that the Alamouti-based protocol enjoys benefits of not only optimal average rate, 
excellent error performance, but also easy decoding in each ARQ rounds. Similarly, one can 
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ask for the LDC to not only satisfy the design criteria presented in this work, but also has 
the "fast-decodable" property [61] for X^") in each ARQ round n. One might have noticed 
that for the LDC-based protocol to work, the channel has to remain constant over different 
ARQ rounds, which is the assumption of this paper. For time-varying channels, protocols 
exploiting both time and spatial diversity should be considered, and how to efficiently design 
the corresponding protocol could be another interesting problem. 



References 

[1] S. Lin, D. Costello, and M. Miller, "Automatic-repeat-request error-control schemes," 
IEEE Communications Magazine, vol. 22, no. 12, pp. 5-17, Dec. 1984. 

[2] D. J. Costello, J. Hagenauer, H. Imai, and S. B. Wicker, "Applications of error-control 
coding," IEEE Trans. Info. Theory, vol. 44, no. 6, pp. 2531-2560, Oct. 1998. 

[3] G. Caire and D. Tuninetti, "The throughput of Hybrid-ARQ protocols for the Gaussian 
collision channel," IEEE Trans. Info. Theory, vol. 47, no. 5, pp. 1971-1988, July 2001. 

[4] S. Sesia, G. Caire, and G. Vivier, "Incremental redundancy Hybrid ARQ schemes based 
on low-density parity-check codes," IEEE Trans. Commun., vol. 52, no. 8, pp. 1311- 
1321, Aug. 2004. 

[5] C. Shen, T. Liu, and M. Fitz, "Aggressive transmission with ARQ in quasi-static fading 
channels," in Proceedings of IEEE International Conference on Communications, May 
2008, pp. 1092-1097. 

[6] , "On the average rate performance of Hybrid-ARQ in quasi-static fading 

channels," accepted to IEEE Trans. Commun., Apr. 2009. [Online]. Available: 
[http:/ /www.ee.ucla.edu/''"congshen/ARQ_TCOM_revise2.pdf, 

[7] L. Zheng and D. Tse, "Diversity and multiplexing: a fundamental tradeoff in multiple- 
antenna channels," IEEE Trans. Info. Theory, vol. 49, no. 5, pp. 1073-1096, May 2003. 

[8] H. El Gamal, G. Caire, and M. O. Damen, "The MIMO ARQ Channel: diversity- 
multiplexing-delay tradeoff," IEEE Trans. Info. Theory, vol. 52, no. 8, pp. 3601-3621, 
Aug. 2006. 

[9] T. Kim and M. Skoglund, "Diversity-Multiplexing tradeoff in MIMO channels with 
partial CSIT," IEEE Trans. Info. Theory, vol. 53, no. 8, pp. 2743-2759, Aug. 2007. 

[10] A. Chuang, A. G. i Fabregas, L. K. Rasmussen, and I. B. Collings, "Optimal throughput- 
diversity-delay tradeoff in MIMO ARQ block-fading channels," IEEE Trans. Info. The- 
ory, vol. 54, no. 9, pp. 3968-3986, Sept. 2008. 

[11] E. N. Onggosanusi, A. G. Dabak, Y. Hui, and G. Jeong, "Hybrid ARQ transmission 
and combining for MIMO systems," in Proceedings of IEEE International Conference 
on Communications, vol. 5, May 2003, pp. 3205-3209. 



38 



[12] M.-K. Oh, Y.-H. Kwon, and D.-J. Park, "Efficient Hybrid ARQ with space-time coding 
and low-complexity decoding," in Proceedings of IEEE International Conference on 
Acoustics, Speech, and Signal Processing, vol. 4, May 2004, pp. 589-592. 

[13] E. Jang, J. Lee, H. Lou, and J. Cioffi, "Optimal combining schemes for MIMO systems 
with Hybrid ARQ," in Proceedings of IEEE International Symposium on Information 
Theory, June 2007, pp. 2286-2290. 

[14] H. Zheng, A. Lozano, and M. Haleem, "Multiple ARQ processes for MIMO systems," in 
Proceedings of the 13th IEEE International Symposium on Personal, Indoor and Mobile 
Radio Communications, vol. 13, Sept. 2002, pp. 1023-1026. 

[15] T. Koike, H. Murata, and S. Yoshida, "Hybrid ARQ scheme suitable for coded MIMO 
transmission," in Proceedings of IEEE International Conference on Communications, 
vol. 5, June 2004, pp. 2919-2923. 

[16] S. Moon, H. Park, A. Goldsmith, and M. Oh, "Bit rearrangement for MIMO retrans- 
missions," in Proceedings of IEEE Clobal Telecommunications Conference, Nov. 2007, 
pp. 3509-3513. 

[17] E. de Carvalho and P. Popovski, "Strategies for ARQ in 2 x 2 MIMO systems," IEEE 
Commun. Letter, vol. 12, no. 6, pp. 441-443, June 2008. 

[18] A. V. Nguyen and M. A. Ingram, "Hybrid ARQ protocols using space-time codes," in 
Proceedings of IEEE 54th Vehicular Technology Conference, vol. 4, 2001, pp. 2364-2368. 

[19] K. Zheng, H. Long, L. Wang, and W. Wang, "Linear space-time precoder with Hybrid 
ARQ transmission," in Proceedings of IEEE Clobal Telecommunications Conference, 
Nov. 2007, pp. 3543-3547. 

[20] S. Wicker, Error Control Systems for Digital Communication and Storage. Englewood 
Cliffs, NJ: Prentice-Hall, 1995. 

[21] D. Chase, "Code combining-a maximum-likelihood decoding approach for combining an 
arbitrary number of noisy packets," IEEE Trans. Commun., vol. 33, no. 5, pp. 385-393, 
May 1985. 

[22] M. Zorzi and R. R. Rao, "On the use of renewal theory in the analysis of ARQ protocols," 
IEEE Trans. Commun., vol. 44, no. 9, pp. 1077-1081, Sept. 1996. 

[23] P. Wu and N. Jindal, "Performance of Hybrid-ARQ in block-fading channels: A fixed 
outage probability analysis," submitted to IEEE Trans. Commun., Nov. 2008. [Onhne]. 



Available: http://arxiv.org/abs/0811.4191 



[24] Y. Liu, V. K. N. Lau, O. Y. Takeshita, and M. P. Fitz, "Optimal rate allocation for 
superposition coding in quasi-static fading channels," in Proceedings of IEEE Interna- 
tional Symposium on Information Theory, Jun.-Jul. 2002, p. 111. 



39 



[25] S. Shamai and A. Steiner, "A broadcast approach for a single-user slowly fading MIMO 
channel," IEEE Trans. Info. Theory, vol. 49, no. 10, pp. 2617-2635, Oct. 2003. 

[26] M. Effros, A. Goldsmith, and Y. Liang, "Capacity definitions of general channels with 
receiver side information," in Proceedings of IEEE International Symposium on Infor- 
mation Theory, Jun. 2007, pp. 921-925. 

[27] E. Telatar, "Capacity of multi-antenna Gaussian channels," AT&T-Bell Labs, Tech. 
Rep., June 1995. 

[28] M. K. Simon, Probability Distributions Involving Gaussian Random Variables: A Hand- 
book for Engineers and Scientists. Kluwer Academic Publishers, 2002. 

[29] J. Hagenauer, "Rate-compatible punctured convolutional codes (RCPC codes) and their 
apphcations," IEEE Trans. Commun., vol. 36, no. 4, pp. 389-400, Apr. 1988. 

[30] J. Ha, J. Kim, and S. W. McLaughlin, "Rate-compatible puncturing of low-density 
parity-check codes," IEEE Trans. Info. Theory, vol. 50, no. 11, pp. 2824-2836, Nov. 
2004. 

[31] N. Varnica, E. Soljanin, and P. Whiting, "LDPC code ensembles for incremental redun- 
dancy Hybrid ARQ," in Proceedings of IEEE International Symposium on Information 
Theory, Sept. 2005, pp. 995-999. 

[32] J. Ha, J. Kim, D. Kline, and S. W. McLaughlin, "Rate-compatible punctured low- 
density parity-check codes with short block lengths," IEEE Trans. Info. Theory, vol. 52, 
no. 2, pp. 728-738, Feb. 2006. 

[33] J. Kim, W. Hur, A. Ramamoorthy, and S. W. McLaughlin, "Design of rate-compatible 
irregular LDPC codes for incremental redundancy Hybrid ARQ systems," in Proceedings 
of IEEE International Symposium on Information Theory, Jul. 2006, pp. 1139-1143. 

[34] C.-H. Hsu and A. Anastasopoulos, "Capacity achieving LDPC codes through punctur- 
ing," IEEE Trans. Info. Theory, vol. 54, no. 10, pp. 4698-4706, Oct. 2008. 

[35] M. Luby, "LT codes," in Proceedings of the 43rd Annual IEEE Symposium on Founda- 
tions of Computer Science, 2002, pp. 271-282. 

[36] A. Shokrollahi, "Raptor codes," IEEE Trans. Info. Theory, vol. 52, no. 6, pp. 2551-2567, 
Jun. 2006. 

[37] E. Soljanin, N. Varnica, and P. Whiting, "Punctured vs ratclcss codes for Hybrid ARQ," 
in Proceedings of IEEE Information Theory Workshop, Punta del Este, Uruguay, Mar. 
2006. 

[38] B. Hassibi and B. M. Hochwald, "High-rate codes that arc linear in space and time," 
IEEE Trans. Info. Theory, vol. 48, no. 7, pp. 1804-1824, Jul. 2002. 



40 



[39] J. K. Zhang, J. Liu, and K. M. Wong, "Trace-orthogonal full diversity cyclotomic space- 
time codes," in Space- Time Processing for MIMO Communications, A. Gershman and 
N. Sidiropoulos, Eds. John Wiley & Sons, 2005, ch. 5, pp. 169-208. 

[40] R. Heath and A. Paulraj, "Linear dispersion codes for MIMO systems based on frame 
theory," IEEE Trans. Signal Processing, vol. 50, no. 10, pp. 2429-2441, Oct. 2002. 

[41] S. M. Alamouti, "A simple transmit diversity technique for wireless communications," 
IEEE J. Select. Areas Commun., vol. 16, no. 8, pp. 1451-1458, Oct. 1998. 

[42] A. Dammann and S. Kaiser, "Performance of low complex antenna diversity techniques 
for mobile OFDM systems," in Proceedings of the 3rd International Workshop on Multi- 
Carrier Spread- Spectrum & Related Topics (MC-SS 2001), Oberpfaffenhofen, Germany, 
Sept. 2001, pp. 53-64. 

[43] J.-C. Belfiore, G. Rekaya, and E. Viterbo, "The Golden Code: A 2x2 full-rate space- 
time code with non- vanishing determinants," IEEE Trans. Info. Theory, vol. 51, no. 4, 
pp. 1432-1436, Apr. 2005. 

[44] M. O. Damen, A. Tewfik, and J.-C. Belfiore, "A construction of a space-time code based 
on number theory," IEEE Trans. Info. Theory, vol. 48, no. 3, pp. 753-760, Mar. 2002. 

[45] A. R. Galderbank and A. F. Naguib, "Introduction to space-time codes," in Space- 
Time Wireless Systems: From Array Processing to MIMO Communications. H. Bolcskei, 
D. Gesbert, C. B. Papadias, and A.-J. van der Veen, Eds. Cambridge University Press, 
2006, ch. 7, pp. 133-153. 

[46] S. N. Diggavi, A. R. Galderbank, S. Dusad, and N. Al-Dhahir, "Diversity embedded 
spacectime codes," IEEE Trans. Signal Processing, vol. 54, no. 1, pp. 33-50, Jan. 2008. 

[47] O. Tirkkonen, A. Boariu, and A. Hottinen, "Minimal non-orthogonality rate 1 space- 
time block code for 3-|- tx antennas," in Proceedings of the Sixth IEEE International 
Symposium on Spread Spectrum Techniques and Applications (ISSSTA 2000), Sept. 
2000, pp. 429-432. 

[48] A. Hottinen and O. Tirkkonen, "Non-orthogonal space-time block code with symbol 
rate two," in Proceedings of Conference on Information Sciences and Systems ( CISS), 
Baltimore, MD, USA, Mar. 2002. 

[49] H. Jafarkhani, "A quasi-orthogonal space-time block code," IEEE Trans. Commun., 
vol. 49, no. 1, pp. 1-4, Jan. 2001. 

[50] C. B. Papadias and G. J. Foschini, "Capacity-approaching space-time codes for systems 
employing four transmitter antennas," IEEE Trans. Info. Theory, vol. 49, no. 3, pp. 
726-732, Mar. 2003. 

[51] R. Ran, J. Hou, and M. Lee, "Triangular non-orthogonal space-time block code," in 
Proceedings of the 57th IEEE Vehicular Technology Conference, Apr. 2003, pp. 292- 
295. 



41 



Texas Instruments, "Double-STTD scheme for HSDPA systems with four transmit an- 
tennas: Link level simulation results," TSG-R WGl document TSGR1#20(01)0458, 
May 2001. 

R. W. Heath and A. J. Paulraj, "Switching between diversity and multiplexing in MIMO 
systems," IEEE Trans. Commun., vol. 53, no. 6, pp. 962-968, June 2005. 

J.-C. Guey, M. R Fitz, M. R. Bell, and W.-Y. Kuo, "Signal design for transmitter 
diversity wireless communication systems over Rayleigh fading channels," in 1996 IEEE 
Vehicular Technology Conference, Atlanta, GA, 1996, pp. 136-140. 

, "Signal design for transmitter diversity wireless communication systems over 



Rayleigh fading channels," IEEE Trans. Commun., vol. 47, pp. 527-537, April 1999. 

V. Tarokh, N. Seshadri, and A. Calderbank, "Space-time codes for high data rate wire- 
less communication: Performance criterion and code construction," IEEE Trans. Info. 
Theory, vol. 44, pp. 744-765, March 1998. 

S. Siwamogsatham, M. P. Fitz, and J. Grimm, "A new view of performance analysis of 
transmit diversity schemes in correlated Rayleigh fading," IEEE Trans. Info. Theory, 
vol. 48, no. 4, pp. 950-956, Apr. 2002. 

J. L. Massey, "Turbo codes, space-time codes, and what next?" in IEEE International 
Telecommunications Symposium (ITS2002), Sept. 2002. 

M. P. Fitz, Fundamentals of Communications Systems. McGraw-Hill, 2007. 

R. A. Horn and C. R. Johnson, Matrix Analysis. Cambridge University Press, 1985. 

E. Biglieri, Y. Hong, and E. Viterbo, "On fast-decodable space-time block codes," IEEE 
Trans. Info. Theory, vol. 55, no. 2, pp. 524-530, Feb. 2009. 



42 



